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rn ■ ABSTRACT 

The RAVE survey, combined with proper motions and distance estimates, can be used 
to study in detail stehar kinematics in the extended solar neighbourhood (solar sub- 
urb). Using the red clump, we examine the mean velocity components in 3D between 
, 6 < i? < lOkpc and —2 < Z < 2kpc, concentrating on North-South differences. 

H ■ Simple parametric fits to the {R, Z) trends for and the velocity dispersions are 

. . .' presented. We confirm the recently discovered gradient in mean Galactocentric radial 

velocity, Vr, finding that the gradient is more marked below the plane, with a Z gradi- 
ent also present. The vertical velocity, Vz, also shows clear structure, with indications 
of a rarefaction-compression pattern, suggestive of wave-like behaviour. We perform 
a rigorous error analysis, tracing sources of both systematic and random errors. We 
confirm the North-South differences in Vr and Vz along the line-of-sight, with the Vr 
estimated independent of the proper motions. The complex three-dimensional struc- 
ture of velocity space presents challenges for future modelling of the Galactic disk, 
with the Galactic bar, spiral arms and excitation of wave-like structures all probably 
playing a role. 

Key words: Galaxy: kinematics and dynamics; Galaxy: solar neighbourhood; 
Galaxy: structure 



1 INTRODUCTION 

The more we learn about our Galaxy, the Milky Way, the 
more we find evidence that it is not in a quiescent, sta- 
* E-mail: mary@aip.de ble state. Rather, the emerging picture is of a Galaxy in 
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flux, evolving under the forces of internal and external in- 
teractions. Within the halo, there is debris left by accre- 
tion events, th e most significant being the Sagittarius stream 
jMaiewski et al . 2003 ) associated wit h the Sagittarius dwarf 
galaxy ( Ibata. Gilmore fc Irwiniri994l ) . In the intersection of 
halo and disk there are also i ndications of accretion debris, 
such as the Aquarius Stream l| Williams et al.l 1 2011) as well 
as loc al halo streams such as the Helmi stream ( Helmi et al.l 
Il999h . While within the disk itself there is evidence of large- 
scale stellar over-densiti es: outward from the Sun there is 
the Monoceros Stream jNewberg et al.l |2002| . lYannv et al.l 
l2003h . inward from the Sun there is the Hercules thick 
disk cloud (iLarsen fc Humphrey^ 1996 , Parke r et al.l |2003| . 
IParker et al.l 2004 . Larsen. Humphrevs fc. 200S ). The origin 
of all such structures is not yet fully understood: possible 
agents for the M onoceros ring are tidal debris (see e.g., 
iPenarrubia et aL | [20 05). t he excitation of the disk caused by 
accretion events fe.g.. | Kazantzidis et al.ll2008l) or the Galac- 
tic warp l|Momanv et al.ll2006l ). while Humphrevs et aLll201ll 
conclude that the Hercules cloud is most likely due to a dy- 
namical interaction of the thick disk with the Galaxy's bar. 

Velocity space also exhibits structural complexity. In 
the solar neighbourhood various structures are observed 
in the distribution of stars in the UV plane, which differs 
significan t ly fr om a smooth Schwarzschild distribution 
l|Dehnenl 1 19981 ). These features are likel y created by 
the combined effects of the Galactic bar ( Raboud et al. 
19981 . iDehnenI Il999l . iMinchey et al.l |2010|. iFamaev et al. 
2007 ) and spiral arm s JPe Simone. Wu &: Tremaine 
2004 : iQuillen fc Minchevl |2005|. lAntoia et al.l |2009| ) or 



both (|Chakrabar'tvr2007. ]Quillen et al.l I2OI1I). D i ssolv 



ing open clusters (jSkulian. Cottrell fc Hearnshawl 1 19971 . 
iDe Silva et aLll2007^ can also contribute to structure in the 
UV plane. Finally, velocity streams in the UV plane can be 
explained by t he perturbative effect on the disk by re cent 
merger events (|Minchev et al.ll2009l . fGomez et al.ll2012l ). 

Venturin g beyond the solar neighbourhood into the 
solar suburb, lAntoia et al. 1 12012| used data from the RA- 
dial Velocity Experiment (RAVE) jSteinmetz et al.l I2OO6I ) 
to show how these resonant features could be traced far 
from the Sun. Additionally, ISiebert et al. I I2OIII (hereafter 
Sll), showed based on RAVE data that the components 
of stellar velocities in the direction of the Galactic centre, 
Vr, is systematically non-zero and has a non-zero gradient 
5{Vr)/5R < —3 kms~^/kpc. A similar effec t was seen in 
the a nalysis of 4400 RAVE red clump giants bv lCasetti et al.l 
I2OIII . The cause of this streaming motion has been variously 
ascribed to the bar, spiral arms, and the warp in conjunc- 
tion with a triax ial dark-matter halo , or a combination of all 
three. Recentlv. ISiebert et al.|[2012l have used density- wave 
models to explore the possibility that spiral arms cause the 
radial streaming, and were able to reproduce the gradient 
with a two-dimensional model in which the Sun lies close 
to the 4:1 ultr a-harmonic resonance of two -arm spiral, in 
agree ment with lQuillen fc MinchevlboOSi and lPompeia et al] 
I2OIII . While structure in the UV plane smears out with the 
increase of sample depth (at d > 250 pc), Gomez et al. 
(2012a,b) showed that the stellar energy distribution pre- 
serves structure associated with dynamical instabilities for 
distances around the Sun as large as 3kpc, consistent with 
the SEGUE and RAVE coverage. 

Two major aims of Galactic Archaeology are to con- 



strain the Galactic potential and determine the Galaxy's 
history of accretion and star formation. These constraints 
can then be compared with cosmological theory with respect 
to the mode in which galaxies are being built-up (see, e.g . , 
IWhite fc Re3l 19781 . lAbadi et al.|[20o3 . iGuedes et al.ll201li ). 
To achieve these goals we must disentangle the complex knot 
of stellar and kinematic over-densities, resonances and veloc- 
ity gradients. 

In this paper we examine the kinematics of the red 
clump giants in the three-dimensional volume surveyed by 
RAVE, focusing particularly on differences between the 
northern and s outhern sides of th e plane. We use the 
Galaxia code l|Sharma et al.l I2OIII ) to create predictions 
what the RAVE survey would be expected to see if the 
Galaxy were a steady-state, axisymmetric object. We take 
advantage of the RAVE internal data set, which contains 
a sufficient number of stars for the stellar velocity field to 
be explored in three spatial dimensions within the region 
6 < R < lOkpc and -2 < Z < 2kpc. We explore the ef- 
fects of systematic and measurement errors, which almost 
everywhere dominate Poisson noise. We consider the effects 
of using alternative sources of proper motion. Systematics 
caused by assuming a single red clump magnitude are also 
explored using the Galaxia model. 

We find that the streaming motion in Vr reported by 
Sll differs North and South of the plane: if Vr vanishes for 
the Local Standard of Rest (LSR) , then there is a large net 
outward motion with a shallow gradient with R above the 
plane, while below the plane, the gradient is steeper and val- 
ues of Vr are lower beyond the solar circle. The behaviour 
of Vz is also complex: in the {R, Z) plane there is a ridge 
of higher Vz that is inclined at 40° to the R axis, cutting it 
just beyond the solar circle. Either side of this ridge, the Vz 
values fall away steeply. Physically, this velocity field implies 
alternate rarefaction and compression, as in a sound wave. 
Thus our three-dimensional v elocity field confirm s the re- 
cent one-dimensional results of lWidrow et al.ll2012l in which 
North-South differences in the velocities of SDSS stars sug- 
gested vertical waves in the Galactic disk. We find that V^ 
is on the whole much more regular, displaying a steady in- 
crease in asymmetric drift with increasing \Z\ and decreas- 
ing R, although some indications of perturbations by spiral 
arm features are present. A simple parametric fit for V4, is 
presented. 

We re-examine the detection by Sll of the gradient in 
Vr using only line-of-sight velocities in light of our three- 
dimensional velocity field. We find intricate velocity varia- 
tions that are sensitive functions of R, Z and 

The paper is organised as follows. In Section [2] we dis- 
cuss our overall approach for the kinematic mapping, the 
co-ordinate systems used and initial selection from RAVE. 
In Section |3] we present a detailed investigation into the use 
of the red clump stars as a distance indicator, where we in- 
troduce a Galaxia model of the red clump to model the dis- 
tance systematics. In Section|4]we present our error analysis, 
investigating systematic error sources first before discussing 
measurement errors, Poisson noise and the final cuts to the 
data using the error values. In Section[5]we examine the 3D 
spatial distribution of our data before going on in Section 
[5] to investigate the variation for Vr, V^ and Vz over this 
space. Section [7] gives the variations in velocity dispersions 
with {R, Z) and a simple functional fit to these trends. In 
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Section [8] we use the line-of-sight method to investigate the 
Vr and Vz trends towards and away from the Galactic cen- 
tre. Finally, in Section [9] we draw our conclusions. 



2 RED CLUMP KINEMATICS WITH RAVE 
DATA 

2.1 General approach 

In this paper we investigate 3D variations of velocities, deriv- 
ing velocities from RAVE line-of-sight velocities (Vios), liter- 
ature proper motions and photometric distances, where we 
focus on the red clump for the distan ce determination. While 
studies such as iPasetto et al.l l2012l lbl. distinguish between 
thick- and thin-disk stars, we do not. Distinguishing between 
the thin, thick disk and halo require either kinematic or 
chemical criteria, each of which have their individual issues 
and challenges, which we wish to avoid. Rather, our aim here 
is to describe the overall velocity structure of the solar neigh- 
bourhood in a pure observational/phenomenological sense. 
The interpretation of this based on thin or thick-disk be- 
haviour is then secondary. 

To visualise 6D phase space we must project it. Different 
projections show various trends more or less clearly. Hence, 
we use a variety of methods to examine the data. We bin 
the data in primarily {R, Z) but also in {X, Y, Z) space (see 
definitions of co-ordinate systems below). For the former we 
will do two projections, one with a contour plot showing 
the mean velocity within the circular ring of given (R, Z) 
- the plotted means are averaged within the {R, Z) plane 
with a moving box car. We also show velocities averaged in 
R within slices in Z. Contour plots in the (X, Y) plane of 
the mean velocities of stars that lie within different slices in 
Z, again employing a moving box car to average in X and 
Y. 

2.2 Co-ordinate systems and Galactic parameters 

As we are working at large distances from the solar posi- 
tion, we use cylindrical co-ordinates for the most part, with 
Vr, and Vz defined as positive with increasing R, <j) and 
Z, with the latter towards the North Galactic Pole (NGP). 
We also use a right-handed Cartesian co-ordinate system 
with X increasing towards the Galactic centre, Y in the di- 
rection of rotation and Z again positive towards the NGP. 
The Galactic centre is located at {X,Y,Z) = (Rq, 0, 0). 
Space velocities, UVW , are defined in this system, with U 
positive towards the Galactic centre. 

To aid comparison with the predictions of Galaxia 
(see Section [2.4|) . we use mostly the parameter values that 
were used to make the model for which Galaxia makes 
predictions. Thus, the motion of the Sun with respect to 
the LSR is taken from ISchonrich. Binnev fc DehnenI |2010| . 
namely Uq = 11.1, V© = 12.24, H^q = 7.25 kms'^ The 
LSR is assumed to be on a circular orbit with circular veloc- 
ity Kirc = 226.84 kms"\ Finally, we take Rq = 8kpc. The 
only deviation from Galaxia's values is that they assumed 
that the Sun is located at Z — 4-15 pc, while we assume 
Z = Opc. 

Sll explored the variation of the observed Vr gradient 
on the values of Vdrc and Rq, finding that changing these 



parameters between variously accepted values could reduce 
but not eliminate the observed effect. We do not explore 
this in detail in this paper but note that the trends that we 
observe are similarly affected by the Galactic parameters; 
amplitudes are changed but the qualitative trends remain 
the same. 



2.3 RAVE data 

The wide-field RAVE (RAdial Velocity Experiment) sur- 
vey measures primarily line-of-sight velocities and addition- 
ally stellar parameters, metallicities and abundance ratios 
of stars in th e solar neighbourhood (Stcinmctz et al."200^; 
IZwitter et al] [2008: Sicbcrt ct al. 2011: Bocchc ct al. 2011,). 
rave's input catalog is magnitude limited (8 < / < 13) 
and thus creates a sample with no kinematic biases. To the 
end of 2012 RAVE had collected m ore than 550,000 spe ctra 
with a median error of 1.2 kms"^ (|Siebert et al.ll2011bl ). 

We use the internal release of RAVE from October 2011 
which contains 434,807 RVs and utilizes the revised stellar 
param eter determination (see the DR3 paper, ISiebert et al.l 
l2011bl for details). We applied a series of cuts to the data. 
Firstly, those stars flagged by the iMatiievic et al.ll20lj au- 
tomated spectra classiflcation code as having peculiar spec- 
tra were excluded. This removes most spectroscopic bina- 
ries, chromospherically active and carbon stars, spectra with 
continuum abnormalities and other unusual spectra. Further 
cuts restric ted our sample to sta rs with signal-to-noise ratio 
STN > 20. lTonrv fc Davis|[l979l cross-correlation coefficient 
R > 5, \fj,a, Hs\ < 400masyr~^, efj,a, e^s < 20masyr~^, 
RV < 600 kms~^ (see Section l4.4p and stars whose Spec- 
traQualityFLAG is null. Where there are repeat observa- 
tions, we randomly select one observation for each star. With 
this cleaning, the data set has 293,273 unique stars with stel- 
lar parameters from which to select the red clump. 

In some flelds at j6| < 25° the RAVE selection function 
includes a colour cut J—K > 0.5 with the object of favouring 
giants. We imposed this colour cut throughout this region 
to facilitate the operation of Galaxia. The selection of red 
clump stars is unaffected by this cut. Additional data cuts 
and selections were performed later in the analysis and are 
broadly a) the selection of red-clump giants (Section 13. 3p . 
b) removal of stars with large extinction/reddening (Section 
I4.1.5|) . and removal of data bins with large errors in mea- 
surements of mean velocity (Section 14. 4|) . 

As a sample with alternative distance determinations, 
we also used an internal release with stellar parameters pro- 
duced by the pipeline that was us ed for the third Dat a Re- 
lease (VDR3) and the method of IZwitter et al.ll2008l firom 
September 2011 with 334,409 objects. We cleaned the sam- 
ple as above, leaving us with 301,298 stars. 

A newer analysis of RAVE data is presented in the 4th 
RAVE data release, Kordopatis et al . (in preparation) . This 
applies an updated version of the iKordopatis et al] 1201 ll 
pipeline to RAVE spectra. Selecting red clump stars us- 
ing the stellar parameters from this pipeline does not pro- 
duce signiflcantly different results to those selected using 
the VDR3 pipeline above; the conclusions of this paper are 
unaffected by the pipeline used for the red-clump selection. 
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Version Mk, Source 

A -1.65 Observation: Alvcs 2000, Grochol ski &: Saraiedini 2002 

B —1.54 Observation: Groenewegen 2008 

C -1.64 + 0.0625|Z(kpc)| Theory: SalarTs fc Girardi 2002 compared to the RAVE population 

Table 1. The normalisations used for the red clump Mk^ magnitudes 



2.4 The Galaxia model 

In lWilliams et al.ll201ll w e introduced the use of the Galaxy 
modelling code Galaxia (|Sharma et al.|[2"oilh to investigate 
the statistical significance of the new Aquarius stream found 
with RAVE. In this study we use Galaxia both to provide 
predictions with which to compare our results, and to inves- 
tigate the effects of contamination of the red-clump sample 
by stars that are making their first ascent of the giant branch 
(see Section [4.1.ip . Galaxia enables us to disentangle real 
effects from artifacts of our methodology, and further to un- 
derstand the population we are examining. 

Based on the Besangon Galaxy model, the Galaxia code 
creates a synthetic catalog of stars for a given model of the 
Milky Way. It offers several improvements over the Besangon 
model which increase its utility in modelling a large-scale 
survey like RAVE, the most significant of which is the abil- 
ity to create a continuous distribution across the sky in- 
stead of discrete sample points. The elements of the Galactic 
model are a star formation rate, age- velocity relation, initial 
mass function and density profiles of the Galactic compo- 
nents (thin and thick disk, smooth spheroid, bulge and dark 
halo). The parameters and function al forms of these com po- 
nents are summarized in Table 1 of ISharma et al.ll201ll . 

We follow a sim ilar methodology to that described in 
IWilHams et"ai] l201ll to generate a Galaxia model of the 
RAVE sample, from which we select a red clump sample fol- 
lowing the selection criteria applied to the real data. A full 
catalog was generated over the area specified by < i < 360, 
5 < 2° and 9 < 7 < 13, with no under-sampling. As de- 
scribed in Section 13.21 the RAVE data were then divided 
into three signal-to-noise regimes and a sample drawn from 
the Galaxia model so that, for each regime, the /-band dis- 
tribution was matched to that of the RAVE sample in 5° x 5° 
squares. The Galaxia /-band was generated after correcting 
for extinction. 



3 THE RED CLUMP 

The helium-burning intermediate-age red cl ump has long 
been seen as a promising standard candle jCannon et al.l 
with its ease of identification on the HR dia- 
gram. In recent years there has been a renewed inter- 
est in the red clump for distanc e determination, e.g., 
(jPietrzvnski. Gieren fc Udalskil [ioo^ ) studied the clump in 
the Local Group, down to the metal-poor Fornax dwarf 
galaxy. Here we investigate the use, selection and modelling 
of this population in solar-suburb RAVE data. 

3.1 The Red Clump /("-band magnitude 

The K-hand magnitude of the red clump, while being rela- 
tively unaffected by extinction, has also been shown obser- 



vati onally to be only weakly dependent on metallicit y and 
age jPietrzviiski. Gieren fc Udalskill200 j : lAlvej|200(]| ). so a 
single magnitude is usually assig ned for all stars in the red 
clump. Studies such as Grochol ski fc Saraiedinil |2002| and 
Ivan Helshoecht fc Groenewegenll2007l have shown that there 
is some dependence on metallicity and age which can be 
accou nted for by the theoretical model of ISalaris fc Girardil 
|2002| . 

If models correctly predict the systematic dependence 
of AIk on metallicity and age, this would have implications 
for our study of kinematics in the solar suburb, for with 
increasing distance above the plane the metallicity will de- 
crease and the population age. From lBurnett et al.ll20lil we 
estimate that the population means change from [M/H] ~ 0, 
Age = 4 Gyr at Z = kpc to [M/H] ~ -0.6, Age = 10 Gyr 
at Z = |4| kpc. According to ISalaris fc Girardill200^ . these 
age and metallicity changes will change the /f-band abso- 
lute mag nitude from Mk(RC) = -LG-fl a.t Z = Okpc to 
Mk(RC) = -1.39 at \Z\ = 4kpc. Hence, the distances to 
stars at higher Z will be systematically underestimated by 
~ 10 per cent. 

Furthermore, there is unc ertainty regarding the aver- 
age value of Mk. lAlvej |2000| gives Mk(RC) = -1.61 ± 
0.03 for local red clump giants with Hipparcos distances 
and metallicities between —0.5 ^ [Fe/H] ^ 0.0. In the 
2MASS system the corresponding absolute magnitude is 
A/k,(RC) = -1.65 ± 0.03. Gr ocholski fc Sa raiedini 200l 
derives a similar value of M_fs-(RC) = —1.61 ± 0.04 for 
-0.5 ^ [Fe/H] =^ 0.0, 1.6 Gyr ^ Age 8 Gyr from 2MASS 
data of open clusters. Extending th e sample of clusters, 
Ivan Helshoecht fc Groeneweger] |2007| derived Mk (RC) = 
-1.57 ± 0.05 for -0.5 [F e/H] 0.4, . 3 Gyr Age ^ 
8 Gyr. The n using the new Ivan LeeuwenI [ioOTi Hipparcos 
parallaxes, iGroenewegenI |2008| found that the Hipparcos 
red clump giants now give Mks(RC) = —1.54 ± 0.05 
(Mif(RC) = -1.50 ±0.05) over -0.9 ^ [Fe/H] sC 0.3, with 
a selection bias, whereby accurate K-magnitudes are only 
available for relatively few bright stars, meaning the actual 
value is likely brighter. The newer values hold better agree- 
ment with the theoretical results of lSalaris fc Girardil [2003 . 
who derived an average value for the solar neighbourhood of 
Mks(RC) = —1.58, derived via modelling with a SFR and 
age- metallicity relation. 

Given the disagreement over the K-hand magnitude for 
the clump and possible metallicity /age variations, we inves- 
tigated the use of the three normalizations for Mk,, for our 
derivation of the RAVE red clump distances. Table [T]sum- 



Here K is used to denote iC-band magnitudes in the 
iBessell fc BrettJi 19891 syste m, while Ks is u sed to denote 2MASS 
values. The relations from ICarpentedr2003l were used to convert 
between the two photometric systems. Note though that J — Kg 
is shortened to J — K in Section 13.21 and beyond. 
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marizes the se values, w here the Version A is the sta ndard 
value from lAlved I2OOOI and iGrocholski fc Saraiedin: 



2002 



2008 



while Version B is the new value from iGroeneweee: 
Version C is derived from the theoretical models and at- 
tempts to take into account some of the possible systematics 
for lower metallicity stars, where we use the variation of age 
and metallicity of RAVE stars described above to develop 
a simple linear relation of Mks with Z, where the value for 
each star is derived iteratively. In Section [4?L2] we examine 
the effect of the different normalizations on our results. 

3.2 Predictions from Galaxia 

To establish the best selection method for the red-clump 
stars we first tried to match the observed distribution of 
the selected RC stars with the prediction of Galaxia. The 
errors in RAVE stellar parameters decrease with STN, so 
the red clump is more localized at higher STN values. Fur- 
ther, the errors in T^s and logg are correlated. To account 
for this in the modelling of the clump, we split the data 
into three regimes: STN values between 20 < STN < 40, 
40 < STN < 60 and 60 < STN. We then generated Galaxia 
models for each STN regime for 5° x 5° squares, match- 
ing the I-band distribution of Galaxia to that of RAVE 
in bins of 0.2 mag. However, as stated in the DR2 and 
DR3 releases, the DENIS I-band photometry has large er- 
rors for a significant fraction of the stars. We therefore used 
Equation 24 from the DR2 paper to calculate I magnitudes 
for those RAVE stars which do not satisfy the condition 
—0.2 < (/denis — J2mass) — (J2MASS — J4^2MASs) < 0.6, as 
well as those that have 2MASS photometry but no DENIS 
values. For the models the reddening at infinity is matched 
to that of the value in Schlegel map and to convert E{B — V) 
to extinction in different photometric bands we used the con- 
version factors in Table 6 of Schlcgcl ct al. 199^. 

Errors were then added to the values of Teff, \ogg and 
J — K from Galaxia to match the distributions in the RAVE 
data. The same seed was used for the random number gen- 
erator for the temperature and gravity values to account for 
the error covariance. Table 2 gives the standard deviations 
of the errors that were added to the values from Galaxia. 
The average J — K error, 0.01, is smaller than the expected 
average observational error in J - A', J2 + eKl = 0.03. 
An error o f 0.01 in J — K correspo nds to an error in TeS 
of ~ 40 7^ (|Alonso. Arribas fc Martinez- RogeiHigQeP . Given 
this small value and that, unlike T^s, the spread m J — K 
does not change with STN, we opted to select the RAVE 
red clump in (J — K, log g). 

3.3 Selection of RC stars 

Figure [1] plots the (J — K,\ogg) plane for the RAVE gi- 
ants within 8 < J < 13 in the three STN regimes along with 
the corresponding Galaxia predictions including errors. Here 
we see the clump aX J — K ^ 0.65 and logg ~ 2.2. There 
are some notable differences between the models and the 
observed distributions. First, the position of the clump in 
log (7 decreases with decreasing STN for the RAVE stars. 
This change is not matched by Galaxia: the predicted dis- 
tribution has mean logg ~ 2.4. With increasing dispersion 
at lower STN it elongates but does not shift. Thus, the ef- 
fect seen in the RAVE data appears not to be astrophysical 



STN regime 


o"(log9) 




a(J - K) 


60 s: STN 


0.25 


25 


0.01 


40 ^ STN < 60 


0.35 


110 


0.01 


20 s; STN < 40 


0.45 


150 


0.01 



Table 2. Standard deviations of the error distributions added to 
the Galaxia model to match the RAVE distributions in the three 
STN regimes. The same seed was used for the logg and T^g to 
mimic the error covariance. 



in origin but rather a result of the parameter determination 
from RAVE spectra (for a discussion on systematic trends of 
stellar parameter estimates with SNR see DR3). Also, there 
is an absence of horizontal-branch stars below J — K < 0.55 
in the RAVE data relative to the prediction by Galaxia. This 
latter difference does not affect our results as they are not 
included in the selection region. 

We select the red clump as those stars within 0.55 ^ 
J - K 0.8, 1.8 ^ logg ^ 3.0. We do not include a STN 
dependence as extending the selection box further up in 
gravity will mean an increased fraction of sub-giants in the 
sample and thus erroneous distances. The red box in Fig- 
ure [1] displays the selection criteria, where we see that for 
20 < STN < 40 some of the clump is cutoff at higher values 
of logg. Applied to our entire cleaned RAVE set, this selec- 
tion yields 78,019 stars. A sample of red clump stars were 
also selected from the Galaxia model in the same way, with 
a slightly smaller number of 73,594 'stars'. 

Figure [5] gives the distributions of the entire RAVE 
data set and red clump stars in /, Vios and overall proper 
motion, /i = \/ + mI- The distributions from Galaxia 
are also given, with an additional spread of 2 kms^^ and 
2.7masyr~^ added to the Galaxia line-of-sight velocity and 
proper motion, respectively, to simulate the observational 
uncertainties. These values give the average value of the er- 
ror in these quantities for the RAVE catalog stars. Here we 
see that Galaxia reproduces quite well the basic distribu- 
tions of the two observational kinematic parameters which 
affirms our use of it in comparisons. The small difference in 
the proper motion distribution could suggest that the actual 
errors in proper motion are slightly larger than estimated. 



4 ERROR ANALYSIS 

The most problematic errors are systematic ones. There are 
multiple ways that these can be introduced into the kine- 
matics. We examine these in some detail before discussing 
measurement errors in Section [4.21 Poisson noise in Section 
14.31 and the final data cuts in Section [4.41 



4.1 Systematic error sources 

4-1.1 Contamination by first- as cent giants 

In the (J — K, log g) plane the red clump is overlaid on the 
first-ascent giant branch. Our selection criteria therefore also 
select first-ascent giants, including some sub-giants, which 
can lead to distance and velocity errors. To assess how this 
contamination affects our distances and kinematics, we use 
the red-clump selected stars from the Galaxia model and 
compute a distance to them assuming a single red-clump 
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Figure 1. The J — K-logg plane for RAVE giants for three STN regimes (a - c) and the corresponding Galaxia models (e-f). Panels 
(d) and (h) show the aggregate RAVE and Galaxia samples respectively. The red clump is visible as an over-density centred at 
{J — K, logg) = (0.65, 2.2) and the selection for the red clump region is outlined by the red box. There is a larger movement of the red 
clump position in log g with decreasing STN for RAVE stars compared to Galaxia. 
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Figure 2. The distributions for I, Vios a^nd overall proper motion ^ for the entire sample of unique, cleaned RAVE stars between 
8 < / < 13 (top) and those selected as part of red clump selected stars (bottom). The / magnitudes for the RAVE stars are from DENIS, 
apart from those with problematic values as discussed in Section 13.21 The Galaxia model reproduces the two observational kinematic 
distributions. 
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magnitude. These results can then be compared to the true 
values for the model stars. Galaxia uses mainly Padova 
isochrones and with a median red-clump K-band magnitude 
of Mk,(RC) = -1.50. This is shghtly different to that ap- 
plied above to our data, but is not significant as we are only 
performing an internal comparison with Galaxia. 

Figure [3^ shows the histogram of the parallax errors. 
The distributions can be decomposed into two Gaussians; 
one representing the red clump and the other first-ascent 
giants. Note that the histogram of distance errors does not 
have such a tidy Gaussian decomposition, hence why we 
work with the parallax errors. The spread in the red-clump 
values is caused by the clump's distribution in age and 
metallicity; the average Mk value chosen is not true for all 
red-clump stars. The first-ascent giants have larger errors 
because while their log g values may overlap with the clump, 
their absolute K-hand magnitude can be quite different. 
From this decomposition we estimate that the ~ 40 per cent 
of selected stars are actually red-clump stars and the rest 
are first-ascent giants. Fortunately however, the mean of 
the background distribution is very similar to that of the 
red clump, with a systematic shift of only ~ 2 per cent in 
distance. In Figure [Sja we plot the cumulative distribution 
of the distance errors. From this we can see that despite 
the high level of contamination, 80 per cent of stars have 
distances errors of less that 25 per cent. 

To establish just how these distance errors affect our re- 
sults, we calculated {Vr, V^, Vz) for stars in the solar cylin- 
der (7.5 < _R/kpc < 8.5) from the pseudo-data produced by 
Galaxia using both the true distance to each star, and the 
distance one infers from the star's apparent magnitude and 
a single RC Mk- The upper panels of Figure 0] show the 
resulting plots of mean velocity components as a function 
of \Z\. In the panels for Vr and Vz the differences between 
the velocities from true distances (solid line) and from ones 
derived from a single absolute magnitude (dashed line) are 
noticeable only at the limits of the surveyed region. In the 
case of the single value of Mk leads to a consistent un- 
derestimation by ~ 5 kms~^. The lower panels of Figure 
|4] show the corresponding results for stars that lie within 
0.5 kpc of the plane, binned in R. Again the impact on Vr 
and Vz of using a single absolute magnitude is evident only 
at the survey limits. We conclude that the use of a single 
value of Mk does not significantly compromise our results. 



^.1.2 Mk normalisation 

In Section [3.11 we introduced three different normalizations 
for the Mk value for the red clump. Figure[5]shows the effect 
of using these different normalizations on the average values 
of Vr, V^, Vz for stars in the solar cylinder (top row of 
panels) and for stars within 0.5 kpc of the plane (lower row). 
We see that the differences in distance produced by the three 
normalizations has an insignificant impact on the average 
velocities: only at the limits of the survey, where sample 
sizes are small, do the differences reach 15 kms~^. Thus the 
dispersion in the velocity errors, on the order of 10 kms~^, 
caused by the single Mk assumption are averaged over and 
thus all but vanish. It is only those bins - the farthest ones 
- with the lower number of stars that show any effect. 

Given that the effect of the absolute-magnitude nor- 



malisation is minor, we use only normalization A for the 
remainder of this analysis. 

^.1.3 Proper motions 

RAVE proper motions are derived from several sources, 
although the maj ority of the catalog gives PPMX 
iRoeser et all 120081 ') or UCAC2 (|Zacharias et al.l |2004| ) 
proper motions. To investigate the contribution of system- 
atic proper-motion errors, we cross-matched with two ad- 
ditional proper-motion catalogs, the third US Naval Obser- 
vator y CCD Astrograph Catalog, UCAC3 (|Zacharias et al.l 
|201(J ) , and th e Yale/San Juan So uthern Proper Motion Cat- 
alog, SPM4 (Girar d et ai]|201ll ). The UCAC3 proper mo- 
tions suffer from strong pl ate-dependent syste matic distor- 
tions north of 5 = -20° (|r oeser et al.ll2010l ). and so we 
exclude stars with 5 > —20° in this catalog. 

Figure[n]shows the mean values of Vr, V^ and Vz with R 
and Z for stars below 5 = —20°, — 1 < K < Okpc calculated 
using the three sources of proper motions, where the restric- 
tions in 5 and Y were made to keep as uniform as possible 
the fractions of stars that have entries in each of the three 
catalogs. The mean velocities do change with proper-motion 
source, with the discrepancies largest between RAVE proper 
motions and those from the SPM4 catalog. The divergence 
is largest for Z < —1 kpc for Vr and V^ with differences up 
to 20 kms~^ in V^. Some of these variations may be the 
result of different coverage of the sample volume - which we 
have tried to minimize with the restrictions on 5, Y above. 
Nevertheless, it does indicate that the proper motions are a 
significant source of systematic error. 

It is difficult to say a priori which catalog is closer to the 
truth; a detailed comparison is beyond the scope of this pa- 
per. Given that the systematic differences between proper- 
motion catalogs can cause significant differences in the de- 
rived velocities, we use all three proper motion catalogs in 
our further analysis, concentrating on the two most diver- 
gent proper motions, namely those in the RAVE and SPM4 
catalogs. 



4-1-4 Binary stars 

Two types of spectroscopic binaries can potentially affect 
our results; single-lined spectroscopic binaries (SBl) stars 
and double-lined spectroscopic binaries (SB2) stars. Both 
kinds introduced an additional velocity variation to the sam- 
ple, with the magnitude much larger for SB2 stars. For the 
SB2 stars, there is the additional problem that the RAVE 
processing pipeline assumes that each observed spectrum is 

that of a single star . 

iMatiievic et al.l l201ll uses repeat observations to es- 
timate the fraction of single-line spectroscopic (SBl) bi- 
nary stars in the RAVE sample, finding a lower limit of 

10 15 percent of the sample consists of SBl stars. The 

technique is biased towards shorter period binaries, but can 
be used to gauge the contribution to the measured disper- 
sions from SBl stars. In Matiievic et al. 201 li the mode of 
the distribution in the velocity variation from the SBl stars 
is 6 kms~^. Undetected, long period binary stars will have 
smaller velocity variations. Thus, the small additional veloc- 
ity dispersion on a limited fraction of the red clump stars 
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Figure 3. a) Histogram of parallax error for Galaxia model stars caused by assuming a single red clump magnitude for stars selected 
in the red clump region. The distribution is decomposed into two Gaussians: one associated with the true RC (dashed blue line) and 
one with the first ascent giants (dotted red line). Their sum is given by the solid green line, b) The cumulative distribution of the 
corresponding distance errors. 
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Figure 4. Comparison between Galaxia model results using the true Galaxia model distances (solid black line) and those assuming 
a RC Mk magnitude (dashed blue line). The average Vr, and Vz for stars within 7.5 < R < 8.5 kpc are plotted as a function of Z 
(top), as well as the values for \Z\ < 0.5kpc as a function of R (bottom). 



from SBl sources can be neglected, especially as we are mea- 
suring mean velocities in this p aper. The fraction of SB2 
stars in RAVE is much smaller; iMatiievic et al.ll201(]| esti- 
mates 0.5per cent of the sample are SB2 stars. Many of these 
are flagged however, as well as any other peculiar stars, and 
have been removed fr om the sample. Unres olved SB2 stars 
are fortunately rare (|Matiievic et al.ll2010l ) and these and 



the remaining SB2 stars do not affect our results substan- 
tially. 

4-1.5 Extinction correction 

As in IWilliams et~al] l201ll . extinction is calculated itera- 
tively from the distances using ISchlegel et al. I Il998l dust 
maps and assuming a Galactic dust distribution as in 
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Figure 5. Comparison between the three version of the red clump Mk normahzation. The average Vn, and Vz for RAVE stars within 
7.5 < R < 8.5kpc are plotted as a function of Z (top), as well as the values for \Z\ < 0.5kpc as a function of _R (bottom). Normalization 
A from Table [T] is given by the solid black line, normalization B by the blue dashed line and normalization C by the green dot-dashed 
line. The results are relatively insensitive to the assumption of a single Mk magnitude. 




Figure 6. As in Figure [S] but a comparison between the results using different proper motion sources for RAVE stars with 5 < 
—20°, — 1 < y < Okpc, with RAVE catalog proper motions (black solid line), UCAC3 proper motions (dotted orange line) and SPM4 
proper motions (dash-dottted red line). 
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Figure 7. Example three stars showing the distribution of the 
MC particles used for error propagation. Some stars, whose ve- 
locities are mostly determined by the proper motions, exhibit a 
large covariance between velocities. 



iBeers et al1l2000l . Note that the maps were adjusted using 
the correction of lYasuda. Fukugita fc Schneider] l2007l . The 
extinction reaches relatively large values of Ak > 0.1 for 
|6j < 10°. However, the effect on the velocities is minor; the 
effect of the extinction correction on the total space velocity 
is greater than 5 kms~^ in only 1.6 percent of red clump 
stars. These stars were excluded from further analysis, leav- 
ing 72,365 stars in our red clump sample. 



4.2 Measurement errors 

The non-systematic sources of errors in the kinematics can 
be broadly broken into (i) the contribution of random un- 
certainties in the measurements and (ii) the finite sample 
size for each volume bin. We will examine each of these in 
turn, before discussing any kinematic cuts that were applied 
to the data. 

The Monte Carlo method of error propagation enables 
an easy tracking of error covariances, so we employ it here by 
generating a distribution of 100 test particles around each in- 



put value of distance, proper motion, line-of-sight velocitj0 
and then calculating the resulting points in {Vr, V^, Vz). 
We assume that the proper motions and line-of-sight veloci- 
ties have Gaussian errors with standard deviations given by 
the formal errors for each star. For the distribution in dis- 
tance for each star however, we generate a double-Gaussian 
distribution in the parallax as in Figure [3l which we then 
invert to derive the distance distribution. This takes into 
account the fact that errors in distance are due both to 
the intrinsic width of the RC and the mis-classification of 
other giants. The ratio between the RC and first-ascent gi- 
ants changes with distance, so we split the data at d = 
1 kpc. The double-Gaussian in parallax then has dispersions 
(0.05, 0.26) mas at a ratio of 2 : 1 for d ^ 1 kpc and 4 : 1 
for d > 1 kpc. Figure [7] shows Monte-Carlo samples of the 
resulting velocity distributions for three typical stars. 

4.3 Poisson noise 

In each of the bins that are used to calculate the mean veloc- 
ities, there are a finite number of stars A'^, so Poisson noise 
contributes to the errors in the derived kinematic properties. 
These errors scale as l/\fN. To establish the error contri- 
bution from this source we used Bootstrap case resampling 
with replacement: for each kinematic quantity we derived 
a distribution in the values by randomly resampling from 
the distribution of values in the bin. The variance in each 
value could then be calculated from the resulting distribu- 
tion. Poisson errors can dominate over measurement errors 
in bins at large distances that contain very few stars. 

To ensure that Poisson errors do not dominate our 
plots, we only use bins which have N stars > 50 and only 
include those points that have errors in the mean of less 
than 5 kms~^ (see Section [4. 4p . 

4.4 Data cuts 

There are several ways in which to prune the data to those 
that are deemed more reliable. However, pruning is liable to 
introduce kinematic biases. We investigated the effects on 
the measured velocities of trimming data via a) proper mo- 
tion, b) error in proper motion, c) distance errors, d) mag- 
nitude of total velocity and e) total velocity error, cVtotai- 
In general, we found that as we increase the cut-off point 
there is a steady increase in velocity dispersion and fluctua- 
tions in the average velocity, asymptotically approaching a 
value as the number of stars approaches the full sample. It 
is therefore diflicult to justify cuts that remove a significant 
proportion of stars. We therefore introduced cuts that re- 
move only the outlier values, as given in Table [S] and we do 
not perform a cut on eViotai- 

Since we seek only to follow trends with R, Z, we do 
not distinguish between halo, thick and thin disk stars. The 
cuts in total velocity therefore aim to be inclusive of halo 
stars; the 6 00 kms~^ limit is ~ 3a the halo dispersion of 
213 kms~^ (|Vallenari et al.|[2006l ). An additional cut in V,p is 
also introduced to limit ourselves to stars that have plausible 
rotation velocities. 

For our analysis we will bin the data in physical space. 



The positions of the stars are assumed to have negligible errors. 
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Cut Reason 

\A^R + (^</> - 220)2 + y|) < 600 kms^^ Remove outlier velocities 

|V0 — 220 1 < 600 kms~^ Remove outlier velocities 

Efi^, e^^ < 20 masyr"^ Remove outlier proper motions 

fia, ^^s < 400 masyr"^ Remove high proper motion stars 

e^/d < 1 Zwitter only, remove large distance error stars 

Table 3. Cuts on the data adopted for this data analysis 



calculating the means and dispersion in each bin. A fur- 
ther data cut was performed post-binning. For each bin we 
calculate the mean (Vr), {V^}, (Vz)- For each of these we 
also calculate an error in the mean, given by standard error 
propagation as 

where q = R, <f>, Z and i — 1,..,N, with the number 
of stars in the average. The values ev,,i are given by the 
MC propagation described in Section 14.21 We remove bins 
with large errors, i.e., e(Vg) > 5 kms~^. This affects only 
peripheral points at large distances from the Sun. 



5 SPATIAL DISTRIBUTION 

Figure [H] gives the XYZ distribution for the RAVE red 
clump stars, where we have differentiated between stars 
above and below the plane for clarity. As in Sll, S12, we 
also plot the lo cation of the spiral arms as i nferred from 
the CO maps of lEnglmaier. Pohl and Bissanta i2008. Going 
from outer to inner (left to right), the arms are the Perseus, 
Sagittarius-Carina and Scutum-Centurus arms. 

Another significant near by feature is the fl e rcules 
thick-disk cloud fHTDO jLarsen fc Humphrevd 1 19961 . 
iParker et al] 12001 iParker et all |2003| detected it via star 
counts in the region I = ±(20° — 55)° both above and below 
the plane at latit udes b = ±(25° — 45)°. This is in Qu ad- 
rant 1. Recently, iLarsen. Cabenela fc Humphrevsll201ll re- 
ported that the HTDC starts at (X, Y) = (0.5, 0.5) kpc in 
our co-ordinates, for 0.5 < \Z\ < l.Okpc. From Figure[8]we 
can see that the RAVE sample in the south int ersects with 
this location of the Hercules thick-disk cloud. lJuric et al.l 
|2008| also found the HTDC in the location (X', ¥', Z') = 
(6.25, -2 - 0, 1 - 2) kpc, corresponding to (X, Y, Z) = 
(1.75, — 2, 1 — 2) kpc in our co-ordiantes. This region is not 
covered by the RAVE survey so we cannot see the northern 
component of the HTDC. In the same paper, another stellar 
over-density was found at (X', Y' , Z') = (9.5, 0.5, 1-2) kpc 
((X, y, Z) = (-1.5, -0.5, 1-2) kpc), which is just missed 
by the RAVE volume. 



6 TRENDS IN MEAN VELOCITIES 

Figure [5] displays the trends in mean Vr, V^, Vz as func- 
tions of R for 0.5 kpc thick slices in Z using 0.5 kpc bins in R. 
Results are shown for three choices of proper motions: those 
in the RAVE, UCAC3 and SPM4 catalogs. Also shown are 
the results of the pseudo-data from Galaxia: the blue line is 
obtained using the true distances to pseudo-stars, while the 



green line uses distances inferred from RC magnitudes. Fig- 
ures [11] and [12] display essentially the same results as grey- 
scale plots in the {R, Z) plane. However, to save space we 
show results only for the proper motions in the RAVE and 
SPM4 catalogs (Figure [TT|) and for the Gatoraa pseudo-data 
(Figure [T2|). Finally, Figure fTSl shows the results obtained 
with the proper motions in the RAVE catalog in full 3D - 
Vr, V^, Vz are averaged over boxes of size 500 x 500 x 500 pc. 
The centres are moved by 250 pc in X and Y and spaced by 
500 pc in Z. 

We now discuss the trends in each velocity component. 
6.1 V^ 

The top two panels of Figure [TT1 show that, as we might ex- 
pect, is largest in the plane and decreases fairly symmet- 
rically both above and below the plane. At given values of 
Z, V4, also decreases inwards. These trends are independent 
of the distances used (RC or Zwitter) and of the adopted 
proper-motion catalog, and they are exactly what the Jeans 
equations lead us to expect: decreases as the asymmet- 
ric drift increases and the asym metric drift increases with 
the velocity dispersion (see e.g. iBinnev fc Tremaind 1 1 9981 . 
eq. (4.34)). Consequently, whether we move inwards or away 
from the plane to regions of increased velocity dispersion, 
will decrease. 

Figure [9] shows that adopting a different proper- motion 
catalog does slightly change the predicted values of V4,: away 
from the plane slightly larger values are obtained with the 
SPM4 proper motions than those in the RAVE catalog, 
while the UCAC3 proper motions give intermediate values. 
However, these differences between results from different 
proper-motion catalogs are much smaller than the difference 
between the observations and the predictions of Galaxia: the 
latter predicts a markedly flatter profile of V^ with Z. In 
other words, the data imply that falls away as we move 
away from the plane significantly more rapidly than Galaxia 
predicts. 

In Section |4]TTT] we saw that the assumption of a single 
RC magnitude can lead to an under-estimation of V^. This 
effect is also evident in Figure [5] in that the dotted green 
lines obtained using RC distances lie below the dashed blue 
lines obtained with the true distances. Since even the green 
lines lie above the data for \Z\ > 0.5 kpc, we conclude that 
the use of a single RC magnitude does not explain the offset 
between the data and the predictions of Galaxia and the 
Galaxy's velocity field must differ materially from that input 
into Galaxia. Comparison of the middle panels of Figures [TT] 
and [12] reveals that, compared to the predictions of Galaxia, 
the observations show more clearly the expected tendency 
for V^ to decrease as we move in a t fixed \Z\. 

In their study of the HTDC, IParker et al] |2004| mea- 
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Figure 8. The XYZ distribution of RAVE red clump gia nts, split into samples above ttie plan e (top) and below the plane (bottom). The 
position of the spiral arms inferred from the CO maps of lEnelmaier, Pohl and Bissant j|2008l are plotted (red lines), giving from left to 
right the Perseus, Sagittarius-Carina and Scutum-Centurus arms. The four Galactic quadrants are also delimited and labelled (Q1..Q4) 
in the central XY plot. Note that the diagonal over-density of stars in Quadrant 4 is purely caused by a greater number of observations 
in that region. 



sured the line-of-sight velocities of thick-disk and halo stars, 
finding that in Quadrant 1 the rotation of this body of 
stars lagged the LSR by 80-90 kms~^, while in Quadrant 
4 the corresponding lag was only 20 kms~^. Far from con- 
firming this effect, the centre-left panel of Figure [TS] sug- 
gests if anything the opposite is true: at X ~ l.Skpc and 
— 1<Z< —0.5 kpc, V4, is lower in the Quadrant 4 than 
Quadrant 1. 

The red curves of Figure [T3] show the positions of spiral 
arms. Close to the plane, at —0.5 < Z < 0.5 there is some 
hint that there is an increasing lag in associated with the 
spiral features. Further from the plane, such an association 
is less clear, as one might expect. 

The [R, Z) dependence of V4, is best described by the 
power-law: 



Vs = ai + [ a2 + a^i 



^kpcy 



z 

kpc 



kms 



(2) 



with Table [4] giving the coefficients for the three differ- 
ent proper motion sources. The for m is not dissimila r to 
that given for hig h- met allicity stars in llvezic et al ] |2008l and 
iBond et al.ll2010l . with however an extra component for the 
R dependence. 



6.2 Vr 

In the middle panels of Figure 1111 contours of constant Vr 
are by no means vertical, whether one adopts the proper 
motions in the RAVE or SPM4 catalogs. In fact, with the 
SPM4 proper motions the contours are not far from hori- 
zontal. Thus the radial gradient SVr/SR = —3 kms^^/kpc 
reported by Sll is at the very least just one aspect of a com- 
plex phenomenon. In fact, if the SPM4 proper motions are 
correct, the trend in Vr is essentially that the further stars 
are from the plane, the more they are moving away from the 
Galactic centre. 

The middle panel of Figure [9] shows this situation in 
a different way by showing that the red curves for the 
SPM4 data are horizontal within the errors at all distances 
from the plane, but they move down from Vr ~ 8 kms"'^ 
at Z ~ —1.25 kpc to zero in the plane and back up to 
~ 8 kms~^ at Z ^ 0.75 kpc. In this figure the black curves 
that join the RAVE data-points tell a different story in that 
in the panels for Z <0 they slope firmly downwards to the 
right. In the two panels for Z > the data points from 
the RAVE proper motions are consistent with no trend in 
Vr with R with the exception of the innermost point in 
the panel for Z ~ 0.25 kpc, which carries a large errorbar. 
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Figure 9. The trends in average velocity as a function of {R, Z) position in the Galaxy for RC stars with RAVE catalog proper motions 
(black solid line), UCAC3 proper motions (orange dashed line) and SPM4 proper motions (red dot-dashed line). Galaxia model results 
are also given, using the real Galaxia distances (dark blue) and RC distances (light green). Error bars give the measurement (thick line, 
short hat) and Poisson (thin line, long hat) errors. The grey open circles in the plot give the result of the fit using RAVE catalog 
proper motions in Section l6.1l 
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Figure 10. As in Figure|9]with however giving the Galactocentric radial velocity with Uq = 14 kms~^. The values are shifted down 
with the overall trends unaffected. 
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Figure 11. The trends in average (top), Vb. (middle) and Vz (bottom) as functions of position in the (i?, Z) plane for RC stars 
with RAVE catalog proper motions (left) and SPM4 proper motions (right). The plotted data are box-car averages over 200 pc X 200 pc 
wide boxes in (R, Z) with 100 pc increments in the co-ordinates of the box's centre. 



The present analysis is consistent with the value of SVr/SR 
reported by Sll in that —3 kms~^/kpc is roughly the aver- 
age of the gradient 5Vr/3R ~ — 7 kms~^/kpc given by the 
RAVE proper motions a.t Z < and the vanishing gradient 
at Z > 0. 

If the dominant gradient in Vr is essentially in the 
vertical direction and an even function of Z as the SPM4 
proper motions imply, the suspicion arises that it is an ar- 
tifact generated by the clear, and expected, gradient of the 



same type that we see in V^. The gradient could be then 
seen to be caused by systematics in the proper motions cre- 
ating a cor relation between the measured valu e of Vr and 
V0 (see e.g. ISchonrich. Binnev fc Asplundll2012l ). In Section 
[8] we re-explore the line-of-sight detection of the Vr, which is 
a proper-motion-free approach to observing the radial gradi- 
ent, which however corroborates the existence of a gradient 
and North-South differences. 

ISchonrichll2012l found a larger value of Uq = 14 kms 
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Figure 12. As for Figure [TT1 but for the Galaxia model using the "true" Galaxia distances (left) and the RC-distances (right). 



compared to that used in this study and Sll, and suggest 
that this larger value would reduce the gradient in Vr. In 
Figure fTOl we plot Vr, the Galactocentric radial velocity with 
Uq — 14 kms~^. The qualitative trends are unaffected by 
the higher value of Uq, with the values only shifted down, 
bringing the values in the 0.5 < Z < 1.0 kpc bin closer to 
those predicted by Galaxia. 

In Figure [13] the location of the spiral arms is overlaid 
for reference on the full XYZ plots of Vr trends. There 
is some indication that structure in Vr is coincident with 
spiral arms in the —1.0 < Z < —0.5 kpc region, with Quad- 



rant 3 and Quadrant 4 showing the largest gradient. With 
a sca le height of the thin disk of 300 pc ([Gilmore fcReidI 
'1983*), at these large Z values we would not expect however 
that spiral resonances would play a role. However, the recent 
results of S12 which explained the gradient in terms of such 
resonances suggests that they do exert an influence at large 
distances from the Galactic plane. The diminishment of the 
gradient at positive Z values is also evident in this flgure. 
Further modelling in 3D is required to understand if the 
North-South differences in the gradient can be explained in 
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Figure 13. The trends in average (top), Vn (middle) and Vz (bottom) as a function of {X, Y, Z) position in the Galaxy for RC 
stars only using RAVE catalog proper motions. The position of the spiral arms (red lines) are also plotted. 



terms of resonances due to spiral structure, which one would 
expect to produce effects symmetric with respect to Z. 

Figure [12] shows the corresponding Galaxia model for 
the red clump. Here we see that there is an absence of the 
observed North-South asymmetry, and further, we confirm 
that the asymmetry is not produced by using the red clump 
distance estimate. We have explored the influence of chang- 
ing the distances using the Zwitter distances. It was found 
that these alternative distance calculations also show the 
same overall R, Z variations for Vr, while exhibiting some 
differences in various regions. This indicates that while the 
distances may influence the quantitative results, qualita- 
tively the overall trends remain the same. 



6.3 Vz 

The bottom two panels of Figure [TT] show that both proper- 
motion catalogs yield a similar dependence of Vz on {R, Z). 
We see a ridge of enhanced Vz that slopes across the panel, 
making an angle of approximately 40° with the R axis and 
cutting that axis about 0.5 kpc outside the solar radius. This 
feature is most pronounced in the bottom-left panel, for 
proper motions in the RAVE catalog. If we assume the LSR 
Vz value is zero, and thus the top of the ridge has Vz = 0, 
we could interpret the stars interior to this point showing 
an overall rarefaction: stars below the plane move down- 
wards, while those above the plane, move upwards. Exterior 
to 71 = 8.5 kpc the behaviour is reversed and shows a com- 
pression; below the plane stars move upwards, while above 
the plane stars move downwards. The amplitude of these 
variations is large; up to 20 kms~^ as seen from Figure (9] 
The results for the SPM4 and UCAC3 proper motions 



do not exhibit quite so large an amplitude in the Vz varia- 
tions, but as shown in Figure[TT]for SPM4, exhibit a similar 
overall behaviour, with the line of higher Vz running diag- 
onal across the RZ plane. However, in Figure [9] we see that 
below the plane, the alternative proper motion results do 
not show the large dip in Vz at = 9 kpc, though this is 
mostly due to the smaller sampling area in the XY plane 
of the alternative proper motion sources. Above the plane, 
the behaviour of the results from the various proper motion 
sources is similar, with clear deviation from the Galaxia 
model's predictions in Figure [T^ The Galaxia model itself 
shows a random pattern of high and low Vz of the order of 
5 kms~^. There is some smoothing spatially of these fluc- 
tuations with the use of RC-magnitude derived distances, 
though the observed pattern is not generated with this dis- 
tance method. In general, the proper motions dominate the 
calculated Vz values, and the differences in observed struc- 
ture reflect this. Nonetheless, it is reassuring that the overall 
pattern is recovered in Figure [11] with each proper motion 
catalog source. 



Recently, IWidrow et al ] |2012l found with main-sequence 
SDSS stars that the vertical number-density and Vz pro- 
flies suggest vertical waves in the Galactic disk excited by 
a recent perturbation. A direct comparison with their re- 
sults and ours is not possible as their stars are outside the 
RAVE spatial sampling region. However, our results support 
this proposition: the rarefaction and compression behaviour 
seen in Figure [TT] is indicative of wave- like behaviour. This 
could then be seen as further "ringing" be haviour of the 
disk c aused by a recent a ccretion event, as in Minchev et al.l 
l2009l : iGomez et a"l]|2012l . Indeed, iGomez et a"l]|2013l further 
suggest that such vertical waves may have been excited by 
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the Sagittarius dwarf as it passed through the disk. Their 
simulations find deviations of up to 8 kms^^ in Vz in a 
heavy-Sgr scenario. 

Turning now to Figure [13] we find some suggestion that 
the Vz field shows some alignment with spiral-arms fea- 
ture, even at significant distances from the plane. Below 
the plane, the Sagittarius-Carina arm at X ~ 1.5 kpc is 
aligned with the low Vz velocities seen with the RAVE cat- 
alog proper motions for y < 0. In this representation, the 
SPM4 and UCAC3 proper motion results also show a dip in 
their Vz velocities around this area, though the feature is 
not as pronounced as for the RAVE catalog proper motions. 
The Perseus arm is also aligned with an area of lower Vz 
at X = -1.25, y < 0, Z < -0.5 kpc. Above the plane, 
there appears to be a reversal of this pattern, with Perseus 
associated with higher Vz velocities for y < 0. Across all 
values of Z, the alignment of the Vz features with the spiral 
arms is less clear for Y > 0, which below the plane may be 
explained by the proximity of the Hercules thick disk cloud. 

These results suggest that there are large-scale verti- 
cal movements of stars, with various cohorts at the same Z 
moving in opposing directions. In .Casetti ct al. 2011, it was 
proposed that while the Galactic warp starts at the solar 
radius, the small elevations of the warp at these distances 
(70 - 200 pc for i? = 8 - 10 kpc (|L6pez-Corredoira et all 
|2002[ )) and the large distances of the RAVE stars from the 
plane would mean tha t kinematics of stars would be little 
affected by the warp. ussiel|[200^ . the nearby Sagittarius- 
Carina arm as traced by star-forming complexes is found to 
lie mostly below the plane by 100 — 200pc, which suggests 
that nearby spiral arms can be influenced by some warp- 
ing. Whether the warp is long-lived or a transient feature 
would have implications for the associated velocities. The 
complexity of the vertical velocity distribution suggested by 
our results would tend to point towards transient features: 
they suggest a non-equilibrium state. A multi-mode trav- 
elling wave caused by a recent perturbation event would 
thus be the most likely explanation for the observed veloc- 
ity structure. 



7 DISPERSION TRENDS 

Figure[T4]displays the trends in the dispersion of V^, Vr, Vz 
as a function of Z for 0.5 kpc thick slices in R using 0.5 kpc 
bins in Z. Also plotted are the results based on UCAC3 
and SPM4 proper motions, and the Galaxia model re- 
sults using the "true" model distances plus RC-calculated 
distances. Unlike the mean-value trends, the differences be- 
tween the various proper motion sources is minimal. We see 
however that the Galaxia model results indicate that the 
use of RC-calculated distances can lead to an increase in 
the measured dispersion values, particularly for Vr. Inter- 
estingly however, there is a reasonable agreement between 
the results of Galaxia and that observed with RAVE. 

We note that the dispersions were corrected for errors 
following the same methodology of lCasetti et al.|[201ll . This 
method adjusts a guess of the true dispersion, ay, in velocity 
component V until 

stddev((l/, - V)/aT) = 1 (3) 
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Table 5. Parameters for the fit given by Equation|4]to the veloc- 
ity dispersion trends using the three proper motion sources; the 
RAVE-catalog compiled results, SPM4 and UCAC3. The param- 
eters are valid for the same ranges given in Table |4] 

where cr^ — ay + eV^, with eVi the MC-derived error in V 
for zth star and V — mea.n{Vi) is the mean over the sample. 

The velocity dispersions exhibit linear behaviour in R 
and parabolic behaviour in Z. Thus, a fit to the dispersion 
trends was found to be provided by the simple form 

R Z ^ 

(yv^bi+b2- i-ba- kms""^ (4) 

kpc kpc 

with Table [5] giving the coefficients for the three velocity 
components for the three different proper motion sources. 
The similar coefficients for each velocity component under- 
lines the similarity of the trends between the three sources. 



8 VELOCITY GRADIENT: LINE-OF-SIGHT 
DETECTION 

In Sll the gradient in Vr was first investigated using the 
line-of-sight velocity, Vios, along narrow lines towards the 
centre and anti-centre. The rationale for this approach was 
to observe the Vr gradient independent of the proper mo- 
tions. The complicated 3D structure of the Galactic Vr pro- 
file was not considered so we revisit this approach, keeping 
in mind the 3D trends. 

8.1 Vr 

Figure 3 from Sll gives the projection onto the plane of the 
mean Vios in bins 200 pc wide as a function of d cos I cos b, 
with the latter a proxy for X in their Galactic co-ordinates, 
with Xsii ~ Rq — dcos I cos b. For \l\ < 5 deg and 175 deg < 
I < 185 deg the values of the radial component, Vr is es- 
sentially the same as the Cartesian value —U. As most 
stars in these narrow cones have cos 6/ sin 6 ~ 1.5, U (and 
so —Vr) is dominated by the term Vios cos 6. In the Ap- 
pendix we list the relevant equations which shows how this 
follows. Sllb compared the mean velocities to those ex- 
pected for a thin disk in circular rotation and adding a ra- 
dial gradient of d{VR)/dR = -3,-5 and —10 kms~^/kpc 
(plus a thick disk lagging the LSR by 36 kms~^), finding 
that the results were consistent with a radial gradient of 
9(Vjj)/9i? = -3 to - 5 kms-Vkpc. 
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Figure 14. The trends in velocity dispersion as a function of {R, Z) position in the Galaxy for RC stars with RAVE catalog proper 
motions (black solid line), UCAC3 proper motions (orange dashed line) and SPM4 proper motions (red dot-dashed line). Galaxia model 
results are also given, using the real Galaxia distances (dark blue) and RC distances (light green). Error bars give the measurement 
(thick line, short hat) and Poisson (thin line, long hat) errors. The open circles give the results of the fitting formula using the coefficients 
from the RAVE catalog proper motions fit. 



From Equations 1 and 7 in the Appendix, we see that 
if we are to consider the Vios as representative of U (and 
so —Vr), technically it would be better to first subtract the 
solar motion from the line-of-sight velocity, to obtain V^^^. 
Otherwise the components above and below the plane will 
be shifted from each other, which can be problematic if the 
sample shifts from above-the-plane to below-the-plane, po- 
tentially producing spurious trends. This was omitted in the 
analysis of Figure 3 from Sll, however, it did not affect the 
conclusions as the models were similarly shifted. 

In Figure [15] we re-examine the trends in V\os cos h for 
the RC giants and stars with Zwitter distances, where we 
plot against R rather than dcoslcosb as this is easier to in- 
terpret and there is little difference between the two values. 
We plot the results with and without the LSR correction, 
with the curve being flattened in the latter with the removal 
of the Uq component. We also plot the trends for a thin- 
disk as in Sll, using 100 Monte Carlo realizations. We see 
here that not correcting for the solar motion does not affect 
the conclusions. Note that with the updated data sets the 
Vr streaming motion is slightly less apparent in both the 
RC and Zwitter results than it was in Sll: the model with 
d{VR)/dR — —3 kms~^/kpc is a better fit to the data. Also, 
the differences between the RC and Zwitter trends around 
the solar radius in Figure [TS] is a result of the different geo- 
metric sampling of the two. 

In Figure [16] we examine the Vios cos b results split into 



North and South components. In this we open up the cri- 
teria a bit to within 7deg of the centre and anti-centre 
direction to reduce the effects of Poisson noise and make 
the trends clearer. Here we see that that differences can 
be discerned between the North and South trends, how- 
ever, they are opposite to that seen in the previous sec- 
tion; the gradient is in the North (—5 to — 10 kms~^/kpc) 
is stronger though absent in the South (0 kms~^/kpc)) for 
R < S kpc. Plotting the corresponding results for Vr in Fig- 
ure [T7] we find that the including the proper motion results 
shifts the overall values so that the northern trends are more 
in line with SVz/5R = —10 kms~^/kpc and the southern 
SVz/SR = —3 to — 5 kms^^/kpc. Despite the disparity with 
the actual numbers, the cause of which is discussed below, 
both Vi'as cos b and the Vr nonetheless exhibit differences be- 
tween the north and south trends, opposite to those found 
in Section [6.21 

To understand the reversal of the trends note first that 
the selected sample cuts across a range of Z as we change 
R so these plots combine the R and Z trends seen previ- 
ously. Second, the necessary restriction on I means that we 
are sampling a very narrow beam and thus do not see the 
global patterns, but only those along that beam. In Figure 
1131 Quadrants 3 and 4 show the largest gradient, which we 
miss with the beams. Hence, these plots emphasize the 3D 
nature of the Vr values in the solar neighbourhood: what 
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Figure 15. As in Figure 3 of Sll, projection of mean RAVE 
Vloa on the Galactic plane in distance intervals of 200 pc to- 
wards the Galactic centre (|(| < 5deg) and anti-centre (175 deg < 
I < 185 deg). Panels (a) and (c) show the results for RC stars 
with heliocentric line-of-sight velocity and corrected to the LSR, 
respectively, while (b) and (d) similarly show the Zwitter dis- 
tance results. The solid curves represent a thin disk population 
with a radial velocity gradient of SVn/SR = 0, —3, —5 and 
— 10 kms~^/kpc, going from green to purple. 
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Figure 16. As in the bottom plots of Figure [TSl but for Galactic 
centre {\l\ < 7 deg) and anti-centre (173 deg < I < 187 deg) and 
split into above (a, b) and below (c, d) the plane. 
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Figure 17. Corresponding Vr trends to Figure 1161 looking at 
the Galactic centre < 7 deg) and anti-centre (173 deg < / < 
187 deg) and split into above (a, b) and below (c, d) the plane. 
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Figure 18. Vj^'cosfe for Galactic centre (|Z| < 7deg) and anti- 
centre (173 deg < Z < 187 deg) directions and split into above 
(a, b) and below (c, d) the plane. The solid curves represent a 
thin disk population with a radial velocity gradient of SVz /SR = 
10, 5, 3, 0, —3, —5 and —10 kms^^/kpc, going from red to 
green to purple. 
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Figure 19. Corresponding Vz trends to Figure 1181 looking at 
the Galactic centre < 7deg) and anti-centre (173 deg < / < 
187 deg) directions and split into above (a, b) and below (c, d) 
the plane. 



you measure very much depends on where you look, be it 
north or south, and at different R and Z values. 

8.2 Vz 

From Equation 6 in the Appendix we can see that Vl cos b 
can be used as a proxy for W = Vz if \ cos b\ > | sin 6| . 
This condition is for the most part met in the low latitudes 
sampled along the |/| < 7 deg and 173 deg < / < 187 deg 
cones used above. So in Figure [TS] we examine the trends 
above and below the plane for Vjcosfc in these directions, 
with Figure [19] giving the corresponding trends in the total 
Vz- Note that we use RC and Zwitter distances here to pro- 
vide the abscissa, supplementing the proper motion data. 
We also plot both positive and negative trends in Vz for 
a thin disk model as above, with values from d{Vz)/dR — 
-10,-5,-3,0,3,5,10 kms'Vkpc. 

Both Figures [18] and [19] show that the trends above 
and below the plane are markedly different, with the 
Vz decreasing with R above the plane at a rate of 
d{Vz)/dR ~ — 10 kms~^/kpc according to the V^ cos b plot 
and —5 kms^^/kpc, to the Vz plot. Below the plane, there 
is a positive trend of d{Vz)/dR ~ -1-5 kms~^/kpc in both 
plots. This positive trend in Vz appears to stop just beyond 
the solar circle at 7? ~ 8.5 kpc. In contrast to Vr above, this 
behaviour is consistent with what was found in Section 16.31 

Note that the differences between the observed magni- 
tude of the trends in Vj^g cos h and Vr , plus V^ cos h and Vz 
can be explained by the fact that both Vi'^^ and V^ contain 
components of U and W. For small I, Vi'^a cos b = U' cos^ b + 
W' cos b sin b and VJ,' cos 6 = — C/' sin 6 cos b + W' cos^ b. The 
cross-terms mean that there is some 'leakage' of the Vz 
trends into Vi^s cos b and Vr trends into 14' cos b, which work 



to either diminish or enhance the observed trend in the single 
component. Nonetheless, the fact that there are differences 
between the north and south for the single components at 
all is further evidence of the 3D variations of the velocity 
values. 

The Vr and Vz trends, as seen via the line-of-sight ve- 
locities and proper motions respectively, are unaffected by 
potential systematics in the distances: a change of the dis- 
tance scale would not affect the fact that a trend is seen 
at all. This is particularly so for the differences observed 
between the northern and southern samples. Furthermore, 
both the SPM and UCAC3 proper motions give similar re- 
sults as in Figures [TS] with a stark negative trend above the 
plane and a positive trend below for the region R < 8.5 kpc. 
Thus, neither the distances nor the proper motions intro- 
duce large systematics into this method of detection. 



9 CONCLUSION 

Using RAVE red clump giants we have examined in detail 
the first moments of the velocity components in a large vol- 
ume around the Sun. We find differences betw een the North 
and South in the streaming motion reported in lSiebert et al.l 
I2OIII in Galactocentric radial velocity, Vr. Above the plane, 
there is a large outward flow with a shallow gradient, while 
below the plane, there are lower values outside the solar 
circle associated with a much steeper gradient in R, partic- 
ularly in Quadrants 3 and 4. 

The behaviour of Vz shows a surprising complexity sug- 
gestive of a wave of compression and rarefaction: there is a 
ridge of higher Vz passing at an angle of 40° to the plane, 
intersecting the plane just outside the solar circle. Assuming 
the LSR Vz to be zero, stars interior to the solar circle and 
above the plane are moving upwards, while those below the 
plane, downwards. Exterior to the solar circle, stars both 
above and below are moving on average towards the plane. 
We confirm these differences by examining the transverse 
velocities along narrow cones towards- and away-from the 
Galactic centre. 

V0 is much more regular than the other two compo- 
nents, showing the most qualitative agreement with the 
mock sample created with the Galaxia model for the 
Galaxy and the expected increase of with increasing Z 
and decreasing R. The measured values are on average lower 
than those given by Galaxia, with the model giving a much 
flatter profile with Z than the data. There is some hint of 
an increased lag in V^ near the plane associated with spiral 
arm features. We present a simple parametric fit to the V^ 
dependence on R and Z. 

We also trace the second moments of the velocities as 
a function of R and Z, providing a simple parametric fit to 
these trends as well. 

The red clump is increasingly being used as a standard 
candle for field stars, given its ease of identification and the 
relative insensitivity of the Mk magnitude to age and metal- 
licity. We modelled our selection of RC stars using Galaxia, 
showing a surprisingly high level of contamination by first- 
ascent giants despite a tight selection in the log g-{J — K) 
plane. However, given that the majority of these giants have 
a similar magnitude to the red clump itself, the effect on the 
distances does not render them unusable. It means though 
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that there is further complexity in the metallicity-age mix- 
ture of selected red clump stars: the population is by no 
means homogenous in age and abundance. 

The assumption of a single Mk magnitude for the RC 
and the proper motions are the largest sources of system- 
atic error in our analysis. Indeed, a deeper study into proper 
motion differences is required to establish which of the cat- 
alogs is most trustworthy. Nevertheless, we have established 
that North-South differences do exist in Vr and Vz despite 
these problems. For Vr a line-of-sight detection, which ex- 
cludes the proper motions, shows gradients above and below 
the plane despite the pencil beams in this analysis pointing 
away from the area of the largest gradient in Quadrants 3 
and 4. These results particularly illustrate the 3D nature 
of the velocity field. For Vz, results using the three proper 
motion sources give the same rarefaction-compression be- 
haviour, albeit with some variation in the details. 

The 3D structure in Vr and Vz presents challenges to 
future modelling of the Galactic disk under the influence 
of the bar, spiral features and any other perturbations (be 
they temporally localised or not). It is not intrinsically clear 
indeed if the structure in the two are coupled or arise from 
different physical mechanisms. 
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motion: 

Vj'jjg = Vias + (Uq cos I COS b + Vq sin i cos b + Wq sinfe), (1) 

V/ = V;-|-(-[/0sin/-|-y0CosO, (2) 

Vl = Vj, + {—Uq cos I sin b — Vq sin I sin b + Wq cos b). (3) 

The reason we do the correction this way round (rather than 
simply in (U, V, W)) is because we wish to see how the line-of-sight 
velocities and proper motions are constituted. For example, the 
Sun is moving up towards the Galactic pole. To correct for this in 
the line-of-sight velocity when looking above and below the plane 
this contribution from Wq is subtracted and added respectively. 
Hence, the solar motion must be corrected for before examining 
trends in Viosi especially if averaging over positive and negative 



The {U, V, W) components are given by 

U = cos I cos b — V/ sin I — cos / sin b 
V = Vj^j, sin I cos b + V/ cos / — sin I sin b. 
W = V{^^sinb+V^cosb, 



(4) 
(5) 
(6) 



where U is positive towards the Galactic centre. For small / this 
reduces to 



U = V{^^cosb-V^smb, (7) 

V = V{, (8) 

W = Vi'^^sinb + V^cosb. (9) 

Similar results are obtained near / = 180 deg, obviously with some 
sign changes. To convert to cylindrical co-ordinates we firstly de- 
fine the Cartesian values corrected for the circular velocity at the 
solar radius; 

Vx = U (10) 
Vy = V + V,,o (11) 
Vz = W (12) 

where we use the nominal value o = 220 kms^^ in the majority 
of this paper. The cylindrical components are then 

VR = ax -Rq)Vx +YVy)/R (13) 
= -ax - Rq)Vy - YVx)/R (14) 
Vz = Vz, (15) 



with R = ^{X - Rq)'^ + y2. For the small i (or i ~ 180 deg) dis- 
cussed in Section [8] where y ~ 0, and for the solar neighbourhood 
where X — Rq < 0, these reduce to Vn ~ —U and ~ V -I- Vc,o 
and so these values are interchangeable in this regime. 
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APPENDIX 

The peculiar motion of the Sun is given in the Cartesian values 
{Uq, Vq, Wq). We define the line-of-sight velocity and proper mo- 
tion vectors in {I, b) co-ordinates corrected for reflex of the solar 



